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Abstract 

The traditional theory of Laplace transformation (TLT) as it was put forward 
by Gustav Doetsch was principally intended to provide an operator calculus for 
ordinary derivable functions of the f-domain. As TLT does not account for the 
behavior of the inverse L-transform at i = its validity is essentially confined 
to t > 0. However, from solutions of linear differential equations (DEs) one can 
discern that the behavior of functions for t < actually is significant. It turns out 
that in TLT several fundamental features of Laplace transformation (LT) are not 
consistently accounted for. To get LT consistent one has to make it consistent with 
the theory of Fourier transformation, and this requires that the behavior of both the 
original function and of the pertinent inverse L-transform has to be accounted for 
in the entire t-domain, i.e., for — oo < t < +oo. When this requirement is observed 
there emerges a new approach to LT which is liberated from TLT's deficiencies and 
reveals certain implications of LT that previously have either passed unnoticed or 
were not taken seriously. The new approach is described; its implications are far- 
reaching and heavily affect, in particular, LT's theorems for derivation/integration 
and the solution of linear DEs. 



1 Introduction 

The mathematical concept which customarily is addressed by the term Laplace Trans- 
formation (LT) was primarily designed as a method for the solution of linear differential 
equations (DEs), i.e., by a kind of operator calculus. Though the theory of LT dates 
back to Leibniz, Euler, Laplace, Petzval, and many 20th-century authors, its presently 
prevalent form, for example in [3 El El, was essentially worked out by Gustav Doetsch 
[H El El SIS]- This version of the theory is in the present article addressed as the traditional 
theory of LT, TLT. 

The theory of Laplace transformation is based on the pair of integral transformations 

L{f{t)} = F{s) = / /(t)e-^* dt; s = a + iuj; a> 0; (1) 
Jo 

L-'{F{s)} = ^{t) = — / F{s)e^' ds, (2) 
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where f{t) is presumed to be a real function of the real variable t. The L-transform F{s) 
is a complex function of the complex variable s = a + iu. In typical applications, t denotes 
time and u denotes circular frequency. 

In TLT the behavior of the inverse transform is customarily characterized by 



at if: = it is left undefined. Indeed, in the realm of ordinary mathematical functions 
- which originally was envisaged by TLT - it is impossible to describe the inverse L- 
transform's behavior at t = 0, except if /(O) = 0; cf. Sect. |21 

By contrast, the original function f{t) is - and has to be - defined at least for < t, 
i.e., including t = 0. If f(t) were at t = insufficiently defined the integral (jT]), i.e., the 
L-transform L{f{t)}, would not exist. 

Hence, in TLT the definition interval of the original function f{t) does not match 
that of the pertinent inverse transform ip{t). In the tradition of Doetsch's theory of LT 
it is widely believed that this mismatch of definition intervals is irrelevant, provided that 
application of Laplace transformation is confined to the interval t > - which in TLT is 
therefore envisaged. However, this assumption is not tenable. For example, the solution 
y{t) of an inhomogeneous linear DE of the type f[y'{t),y"{t), . . .] = x{t), even when 
envisaged only for t > 0, in general depends on the behavior of x{t) for — oo < t [T^ . 
Thus, when the solution y{t) is to be obtained by LT it is crucial which kind of function 
is by L{x{t)} actually represented in the interval — oo < t. This fact is in Sect. lA.ll 
illustrated by an example. 

The most obvious symptoms of TLT's inconsistency may be listed, as follows. 

a) The mismatch of definition intervals appears to disallow concatenation of L-transforms. 
As the inverse transform (p{t) is at t = undefined it can, rigorously, not be allowed to 
be subject to another L-transformation. The fact that concatenation turns out actually 
to be possible does not make TLT formally consistent in this respect. In Sect. 12.41 it is 
explained why concatenation of L-transforms is possible. 

b) The mismatch of definition intervals appears to exclude impulse functions at t = - 
such as 6{t) - from L-transformation. The delta-impulse 6{t) is only at t = different 
from null; thus, according to TLT 6{t) can not exist in the inverse L-transform. The fact 
that the impulse functions actually can be retrieved from their L-transforms does not cure 
this formal inconsistency of TLT. In Sect. 12.21 it is explained why impulse functions are 
LT-consistent. 

c) TLT's derivation theorem 



for t < 0, 
f{t) for t > 0; 



(3) 



L{f{t)} = sL{m}-m 



(4) 



is in conflict with the definition interval of TLT, i.e., t > 0. The real constant /(O) that 
in (0} appears in the L-domain represents a t-domain function of its own, namely, the 
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delta impulse f{0)S(t) at t = 0. This impulse is outside TLT's definition interval and 
therefore should be regarded as irrelevant [12]. Inconsequently, TLT praises the "initial 
value" /(O) as one of its most advantageous features, as that value eventually appears in 
the solutions of linear DEs and is utilized to account for a system's initial state at t = 0. 
By its inconsequent treatment of the initial value TLT implicitly admits that the behavior 
at t = of functions and of their L-transforms actually is significant. 

d) TLT's derivation theorem (j^ is inconsistent with TLT's integration theorem 

L{f{t)} = -L{f{t)}. (5) 

s 

The two theorems differ by the initial value /(O). In TLT, this unexpected and unex- 
plained discrepancy is customarily tolerated. In Sect. 13.41 the relationships between the 
theorems for derivation and integration are outlined and the origin of the conflict between 
Eqs. (jH) and © is revealed. 

e) TLT does not in general keep its promise to provide the solution of linear DEs, i.e., 
for t > 0. There are discrepancies involved which tend to be disguised by formal pseudo- 
consistency. This notion just restates what was said above about solution of linear DEs. 
In Sect. lA.ll an example is described of this kind of failure. 

f) In particular, TLT's general solution of the linear DE suffers from the so-called initial- 
value conflict. In TLT it has become customary to work around this conflict by "patching" 
the original solution. The initial- value conflict is discussed and explained in Sect. 14.41 

These observations indicate that in TLT certain fundamental aspects of LT's behavior 
are not consistently accounted for. The mismatch of definition intervals and its conse- 
quences need to be resolved rather than circumvented and/or ignored. Laplace transfor- 
mation needs to become updated. 

The present article offers an outline of a new, alternative approach to LT. The prob- 
lem of LT-consistency of t-domain functions is re-inspected and resolved. From the results 
such obtained there emerge new methods for the solution by LT of both the linear inho- 
mogeneous and the linear homogeneous DE. 

Crucial results and observations are noted and emphasized as theorems. Several fa- 
miliar theorems of LT, such as, e.g. the fundamental 

superposition theorem 

L{cMt) + C2f2{t) + . . .} = L{ci/i(t)} + L{c2/2(t)} + ... 

= ciL{/i(t)} + C2L{/2(t)} + ..., (6) 

are not affected by the new insights. Another group of theorems become restated and 
re-justified without assuming a new form. A third group includes theorems that become 
more or less drastically modified as compared to their familiar form. Eventually, there is 
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a fourth group, i.e., of theorems which may be regarded as new - at least in so far as in 
TLT they do not play a role. 

Several complementary explanations and examples are exiled into an appendix. This 
article is based on, and complements, earlier related work of the present author P ^ IT ^ ITH]. 



2 Getting Laplace transformation consistent 



TLT's heel of Achilles lies at t = 0. For LT to be consistent it is not sufficient that 
ip{t) = f{t) for t > 0; rather, the condition 

^(0) = /(O) (7) 

must also be fulfilled. To suggest the implications of this requirement, the behavior of 
the inverse L-transform pertinent to a continual function f{t) is illustrated in Fig. 1. 





no) 




^(0) 




\f{t) 


Q. 


\>{f] 


















f 



Fig. 1. Continual function f{t) (left) and the per- 
tinent inverse L-transform ip(t) (right). ip(t) is 
defined for — oo < i and includes at t = the 
connect interval that extends from to /(O) 



From the figure it becomes apparent that for functions which at t = assume a 
definite unique value /(O) ^ the inverse L-transform ip{t) includes at t = an abrupt 
transitional section, i.e., from to /(O). As a consequence, such type of function can 
not meet the criterion ((7j), because a section of (p{t) at t = can not be equal to the 
definite unique value /(O). This kind of non-LT-consistency applies, in particular, to the 
prominent class of continual derivable functions. By contrast, if f{t) is a priori defined 
as a causal function, i.e., f{t) = for t < 0, the criterion may actually be met, i.e., for 
certain conditions which will be discussed below. 

For brevity and simplicity, the discussion of LT-consistency is in the present article 
focussed on the dichotomy between causal functions and d-functions. The term d-function 
denotes the class of bilateral continual functions that are in the ordinary sense derivable 
as many times as required. It should be kept in mind, in particular, that any solution of 
a linear homogeneous DE of finite order is a d-function. 

To achieve an explicit account of LT's behavior at t = 0, the definition interval must 
obviously be expanded from t > into t < 0. As the scope of LT's formula for inverse 
transformation Q already encompasses the entire t-domain, it is the scope of the formula 
for L-transformation (Q) which has to be expanded. To accomplish this kind of expansion 
it is not required to challenge the definitions of LT; the bilateral scope is already implied 
in Eq. (P). 
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2.1 Laplace transformation by Fourier transformation 

The implications of the latter notion become apparent when one exploits the intimate re- 
lationship that exists between Laplace transformation and Fourier transformation. Equa- 
tion is equivalent to the following "unilateral" Fourier-transformation formula for the 
function f(t)exp{—at), i.e., 

POO 

L{f{t)} = F{s) = Fiou, a) = / /(t)e-'^*e--* dt. (8) 

Jo 

However, there is not really such a thing as unilateral Fourier transformation; Fourier 
transformation is inherently bilateral. The actual analysis interval (the "scope") of the 
transformation (jH)), and thus of (PJ), is not determined by the integral's limits but by the 
reciprocal of the frequency spacing df = du;/(27r) of the corresponding Fourier- integral 
representation ^J^j. The actual analysis interval of both (PJ) and (jHl) encompasses the 
entire t-domain — oo < t < +oo, and the low limit t = of the integral ([Q) corresponds to 
the center of the analysis interval. Though these implications of Fourier transformation 
are fairly elementary, apprehension of them appears to be scarce; cf. Sect. IA.2I 

Hence, the unilaterality of the integration interval in both (jH)) and is not equivalent 
to unilaterality of the analysis interval but rather indicates causality of the transformed 
function. Equation (jH} has to be equivalent to ordinary, i.e., bilateral Fourier transforma- 
tion of a causal, i.e., bilaterally defined function fdt) exp(— at), such that 

L{fm = F{u, a) = / /(t)e-*e-i^* dt = / /,(t)e-'^*e-'-* dt. (9) 

Jo J-oo 

For LT to be consistent with Fourier transformation it is necessary that the causal 
function fdt) be defined in such a way that the two integrals in Q are fully equivalent. At 
first sight this requirement is met, e.g., by the definition fdt) = for t < 0; fdt) = f{t) 
for t > 0. However, this definition of the causal function is not actually sufficient, because 
it leaves the transition from /c(— 0) to /c(+0) undefined. There is another condition 
involved: The second (bilateral) integral in Q has to be consistent with the pertinent 
Fourier- integral representation, i.e., the "inverse Fourier transform" 

at <.+oo 

^(t) = L-^{L{f{t)}} = — / F{uj, (T)e^^* dto for - oo < t. (10) 

ZTT J-oo 

As is suggested in (fTUI) . this expression is equivalent to Q. The Fourier- integral 
representation (fTU]) of (p(t) is continuously defined for — oo < t < +oo. In particular, while 
ip{t) is causal, the transition from (p{—0) to v^(+0) is not undefined but either includes 
the so-called connect interval (see below, in particular Sect. 12. 3|) . or its derivative(s), i.e., 
the impulse functions S^"'\t) {n = 0,1,...). The Fourier integral's evident capability 
to represent impulse functions proves that neither the inverse Fourier transform nor the 
inverse L-transform are at t = undefined. As a consequence, for ip{0) = /c(0) to hold 
/c(0) must not be left insufficiently defined. In TLT this requirement is ignored. 
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Below, the consistent definition of fc{t) and tlius tlie consistent expansion of LT's 
transformation formula (Q), (jH)) into an equivalent bilateral form, is approached by three 
steps, advancing from impulse functions to the unit step function and finally to the entire 
class of causal functions. 

2.2 Impulse functions 

The requirements just outlined for fc{t) are a priori met by the delta impulse and its 
derivatives. Thus, Eq. (0) applies to fc{t) = S^^\t) and one obtains 

= s"; n = 0,l,... (11) 

As the impulse functions are only at t = different from null, inverse transformation 
is by (fTUI) achieved for exp(st) = 1, i.e., by inverse Fourier transformation of the function 
(iu;)". One obtains 

L-'{L{S^-\t)}} = S^'^\ty, n = 0,l,... (12) 

It is thus established by the theory of Fourier transformation that both the L-transform 
of 5(")(t) and its inverse L-transform exist, and that the latter is identical to 6^'^\t). The 
impulse functions are LT-consistent. 

As (5'^")(t) is defined for — oo < t < +oo, while (5'^"^(t) = for t 7^ 0, existence 
of the L-transform L{5^"\t)} is consistent with the condition that for the "unilateral" 
transformations (jT]) and (jHI) to apply f(t) must be defined for t > 0. Thus, even when 
one sticks to the "unilateral" form of the transformation there is no conflict. It is only 
the inverse L-transform of impulse functions that is in conflict with TLT's assumptions. 

2.3 The unit step function and the connect function 

Also the unit step function is causal and may be regarded as just another member of the 
class of impulse functions. The unit step function can be defined as the integral of the 
delta impulse at t = 0, i.e., by 

u{t) := f 6{t) dr; - 00 < t < +00. (13) 
J —00 

Equation (fT^ may be regarded as an implicit definition of u(t). When the unit step 
function is explicitly defined it is important to preserve its definition at t = 0. This can 
be accomplished in the form 

u{t) = for t < 0, 

= Uo{t) for t = 0; Uo{t) = {0 ... 1}, 

= 1 for t > 0. (14) 
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This definition accounts for tlie abrupt transition from to 1 that occurs at t = 0, namely, 
by inclusion of the so-called unit connect function, uo(t). The pseudo-function uo{t) may 
be conceived of as an infinite set {0 ... 1} of real numbers (a distribution) that exists at 
t = 0. 

Inclusion of uo{t) in u{t) is indispensable not only for LT-consistency of u{t) but also 
for u{t) to be consistent with the concept of the delta impulse. LT-consistency of u(t) 
requires that the condition ((7j) is met, i.e., that L~^{L{u(t)}} = u{t) at t = 0. When 
the interval extending from «(— 0) = to m(+0) = 1 is left undefined, fulfilment of the 
condition ((Tj) remains undecided. 

The existence of the delta impulse, which is conceptualized as the first (non-ordinary) 
derivative of u{t), crucially depends on existence of the connect function Uoit). With 
respect to (jl4|) there holds 

5{t)=u',{t)=u'{t). (15) 

Thus, indeed, inclusion of Mo(^) in u{t) is indispensable. One can not in earnest conceive 
the delta impulse to be the derivative of a gap. 

As the unit step function is the integral of the delta impulse, it follows from the 
results depicted in Sect. 12.21 that the unit step function as defined by (fT5|) and is 
indeed LT-consistent; i.e., there holds 

L-^{L{u{t)}} = u{t) for - oo < t < +oo. (16) 

Notably, the L-transform of the unit step function is equivocal. There holds 

L{uit)} = L{1} = 1/s. (17) 

With respect to the interval t > the functions f{t) = u{t) and f{t) = 1 differ by the 
connect function Uo{t) which is included in f(t) = u{t) but not in /(t) = 1. The connect 
function does not become explicitly apparent in the L-transform. As a consequence, from 
the L-transform 1/s one can not tell whether it was obtained from f{t) = u{t) or from 

fit) = 1- 

From the fact that both u{t) and its derivatives (the impulse functions) are LT- 
consistent one can conclude that uo{t) invariably is contained in the inverse L-transform. 
According to (fT^ ITTj) there holds 

L-'{L{u{t)}} = L-'{L{l}} = u{t), (18) 

and u{t) contains UQ{t). Thus, for a real function f{t) to be LT-consistent it is indis- 
pensable that it contains the connect function. This is why f{t) = u{t) is LT-consistent 
whereas f{t) = 1 is not. 

The somewhat confusing behavior of 'Uo(^) emerges from the fact that uo{t) is a null- 
function, i.e., 

/ uo(t) dt = 0; T > 0; (19) 

J-T 
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and, therefore, 

L{uo{t)} = 0. (20) 

Whereas uo{t) as such does not become apparent in L-transforms, its derivatives do. 
The derivatives 

=m(")(^) =5^""'^^); n = l,2,... (21) 

have the L-transforms 

L{4"^(t)} = L{M(")(t)} = L{5("-^)(t)} = n = l,2,... (22) 

and these are different from nulL 

The occurrence of the connect function in the inverse L-transform was imphcitly 
pointed out, e.g., by Doetsch [21111 El- He proved that the inverse L-transform is identical 
to the original function except for a null-function, i.e., a function whose integral is null. 
Thus, Doetsch in effect anticipated the involvement of the connect function. However, the 
fact that the null function does not become apparent in L-transforms led him to define 
inverse L-transforms that differ only by a null-function to be identical. As a consequence, 
in TLT the connect function is being ignored. This is a serious mistake, i.e., for the 
following reasons. 

a) Without inclusion of uo(t) one can not obtain a solid conceptual definition of LT- 
consistency, as was pointed out above; cf. Fig. 1. 

b) Without inclusion of uo{t) one can not obtain a solid formal definition of LT-consis- 
tency, because the unit-step redundancy theorem u{t) ■ u^'^\t) = u'^"'\t) does not hold; cf. 
Sects. 1231 

c) Without inclusion of Uo(^) neither the derivatives of u{t), i.e., the impulse functions 
^''"'•'(t), nor their L-transforms are sufficiently defined. 

2.4 Causal functions 

Utilizing the consistent definition (fT^ IT^ of the unit step function, finally the causal type 
of function fc{t) can be defined in the familiar way, i.e., 

f,{t) = u{t)f{t). (23) 

It is inclusion of the unit connect function UQ{t) in the unit step function u{t) that 
makes the expression (j23j) universal. Equation (|23|) holds irrespective of whether or not 
/(t) itself is causal. f{t) may be either a d-function /d(t); a causal function of the form 
[w(t)/d(t)]; or an impulse function 5^'^\t). This follows from the identities outlined in 
Sect. lA. 31 which can be subsumed by the 
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unit-step redundancy theorem 

u{t) ■ u^''\t) = u^''\t)] n = 0,l,... (24) 



By this theorem Eq. (j25|) remains in effect unchanged when both sides are multiplied by 
u(t); in particular, there holds 

uit)m = at). (25) 

Multiplication by u{t) of any kind of causal function is redundant. 

On the basis of these observations the bilateral formula for L-transformation envisaged 
in (0) can be restated, and its equivalence to ((T)) can be expressed by the 



bilaterality theorem 



r+oo 

L{f{t)} = I u{t)f{t)e-'' dt. (26) 



It is the LT-consistent definition (fT^ of u{t) which renders Eq. ()2(i|l consistent, i.e., by 
consistency with Fourier transformation; and Eq. becomes consistent with Fourier 
transformation by its equivalence to (j^Uj) . The mathematical implications of (0) are the 
same as those of 



While the above three-step approach to the formula (|26p is helpful by its elucidating 
implications, it should be noticed that formally the expansion of (Q) into (f^^ can be 
obtained by one single step, namely, 

L{fm= /(t)e-^* dt = hm / r(t,T)/(t)e-^*dt, (27) 

Jo ^ — J —oo 

where r{t,T) denotes the unit ramp function of which an example is illustrated in Fig. 2. 



r{t] 









T 

/ — 



Fig. 2. The unit ramp function r(t) as an approximation to the unit 
step function u{t). For T — > u{t) emerges from r(t), and the unit 
connect function UQ{t) emerges from ro(t) 



By letting T — the unit step function u{t) emerges from the unit ramp function 
r{t,T), and the connect function UQ{t) emerges quite naturally from the ascending part 
ro(t) of the unit ramp function. 

From this approach it becomes apparent that the transition from u{—0)f{t) = to 
u{+0)f{t) = /(+0) is not empty but is a vertically ascending continuous function, i.e., 
uo{t). From this notion there emerges another definition of uo{t), namely, 

Mo(t) = limro(t,T). (28) 
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The function ro(t) does not necessarily have to be hnear - such as in Fig. 2 - but may 
be represented by any kind of real continuous function that in the interval < T rises 
monotonically from to 1. 

Comparison of to reveals a crucial feature of LT which in TLT is ignored, 
namely, LT's inherent ambivalence. There holds the 

ambivalence theorem 



From an L-transform F{s) one can not tell whether it was obtained from f{t) or from 
u{t)f{t). This kind of ambivalence corresponds to the ambivalence of L{u{t)} that was 
noted by (|17p. In both cases the ambivalence is significant by the presence of UQ{t) in u{t) 
and by the fact that L{u[,"^(t)} = s^-^ ^ (n = 1, 2, . . .). 

As a consequence of the equivalence of (f^Bj) to the formula for Fourier-transforma- 
tion, the inverse L-transform obtained by is invariably and unequivocally identical to 
u{t)f(t). There holds the 

causality theorem 

ip{t) = L-^{L{f{t)}} = L-^{L{u{t)f{t)}} = u{t)f{t) for - oo < t < +oo. (30) 
The inverse L-transform is the causal companion of the original function f{t). 

The latter two theorems warrant consistency of concatenation of L-transformations. 
From (123) and ^ there follows the 

concatenation theorem 



In TLT this theorem can not exist because 93(0) is undefined such that, rigorously, L{ip{t)} 
can not be supposed to exist. 

As the behavior of the L-transform L{/()f:)} is in every respect characterized by the 
causal inverse L-transform v'(t), in a sense the L-transform of f{t) actually constitutes 
the L-transform of the causal function u{t)f{t). In many contexts - e.g., that of the 
derivation/integration theorem (cf. Sect. El) - it is indeed quite helpful to observe the 

alias theorem: 

The notation L{/(t)} of an L-transform ultimately is an alias for L{u{t) f {t)} . 

From the above insights there emerges a concise proof of the 
convolution theorem 



L{f{t)} = L{u{t)f{t)} = F{s). 



(29) 



L{ip{t)} = L{f{t)}. 



(31) 



L{/i(t)} ■ L{Mt)} = l{ h{t - t)Mt) dr} 



(32) 
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The proof of (jH^ can be based on the identity 

/ giit)dt-[ g^it) dt = u{t) [ [ g,{t-T)g2{T) dr dt; T >0; t<T, (33) 

J~oo J— CO Jo Jo 

which holds for any pair of causal functions gi^2{t)- Letting gi,2{t) = u{t)fi 2{t) exp{—st) 
and T ^ oo, one obtains from (jHHj) 

rco rco rco rt 

/ /i(t)e-*dt-/ f,{t)e-^Ut= e-^*/ /i(t-r)/2(r)drdt, (34) 
Jo Jo Jo Jo 

where the "factor" u(t) is appropriately accounted for. By ((T)) Eq. (j2D), indeed, is equiv- 
alent to dni. 

2.5 Testing LT-consistency 

For a function f(t) to be LT-consistent it is with regard to (jHUj) crucial that 

n(0)/(0) = /(O). (35) 

In this equation neither ■u(O) nor /(O) necessarily denote unique definite values. The 
symbol m(0) is just a synonym for the pseudo-function Uo(t). As f(t) may be any linear 
combination of d-functions, causal functions, and derivatives, /(O) may actually denote 
any linear combination of definite unique values, connect functions, and derivatives of 
the latter, i.e., impulse functions. Keeping this in mind, and observing the unit-step 
redundancy theorem (j211), the criterion for LT-consistency of any type of real function 
f{t) can be universally expressed by the 

LT-consistency theorem: The function f{t) is LT-consistent if, and only if, 

u{t)f{t) = f{t) for - oo < t < +00. (36) 

The criterion (jHUj) provides for the formal proof of the conclusion that was already 
drawn from Fig. 1, namely, that d-functions /d(i) are non- LT-consistent: u{t)f^{t) ^ fd{t)- 
By contrast, causal functions fc{t) are LT-consistent, because by the unit-step redundancy 
theorem there holds u{t)fc{t) = fc(t). 

A useful application of the LT-consistency theorem is verification of LT-consistency of 
shifted functions. From the formula for inverse L-transformation Q one obtains in the 
familiar way the relationship 

v{t-T) = L-'{L{f{t)}.e-n, (37) 

which indicates that multiplication of L{f{t)} by exp(— sr) shifts the causal inverse trans- 
form by any positive or negative amount r. Utilizing (j30|l one obtains from (j37jl by 
L-transformation the 
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shifting theorem 

L{f{t)}-e~'^ = L{u{t-T)f{t-T)}; - cx)< r < +00. (38) 

This expression holds for any type of real function f{t) and for any r. However, it is only 
for r > that ()38|) is consistent with the alias theorem, indicating that ()38|) is only for 
r > LT-consistent. Indeed, for a shifted causal function fc{t — r) the criterion (j36|) 
reads 

M(t)-/e(t-r) = /c(t-r), (39) 

and this condition is met only for r > 0. 

The shifting theorem makes particularly apparent that Laplace transformation virtu- 
ally is confined to causal functions. TLT's attempt to express the shifting theorem in 
terms of ordinary derivable functions is awkward and incoherent IHEl- 

Yet, Laplace transformation of non-causal functions, in particular, d-functions, does 
not entirely have to be ruled out. L-transformation of d-functions can be sensible and 
useful; cf. Sects. 13. 2| 13. 4| 14.21 One just has to keep in mind that the L-transform of a 
d-function, L{fd(t)}, actually is the L-transform of u(t)fd{t), i.e., L{u(t) fd(t)} . 

2.6 Summary 

Laplace transformation of a real function f{t) is said to be consistent if the pertinent 
inverse L-transform is identical to f{t). For t > LT is in this sense consistent for any 
type of real function that can be L-transformed at all. However, such confined consistency 
is not sufficient for LT to provide a coherent system of operator calculus. Even in the 
realm of TLT the behavior of functions at t = turns out in effect to be involved. Explicit 
inclusion of the point t = into LT's deffnition interval requires 

a) appreciation of the implicit bilaterality of the L-transformation formula (0); and 

b) mathematical description of the behavior at t = of both f{t) and f{t), i.e., by the 
connect function uo{t) and/or its derivatives. 

From a) there emerges a bilateral equivalent of the L-transformation formula that 
makes LT compatible with the theory of Fourier-transformation. This formula in turn 
is dependent on consistent deffnition of the unit step function, i.e., according to b). 
L-transformation turns out to be ambivalent, i.e., there holds L{f{t)} = L{u{t)f(t)}. 
Inverse L-transformation is unequivocal; the inverse transform has the form u{t)f{t). 
Therefore, only causal functions are LT-consistent. 
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3 The derivation/integration theorem 

There are two fundamentally different approaches to obtaining LT-theorems such as those 
for derivation/integration: 

a) Determination by (Q) of the L-domain operation that corresponds to derivation/inte- 
gration of the original function f{t). 

b) Determination by (j21) of the L-domain operation that corresponds to derivation/inte- 
gration of the inverse L-transform f{t). 

L-transformation by ([Q) or is ambivalent, whereas inverse L-transformation by 
(j21) is unequivocal. Moreover, any kind of application of LT, in particular, to the solution 
of linear DEs, is in the first place dependent on the inverse L-transform ip{t) - as opposed 
to the original function f(t). Therefore, only the approach b) is adequate. 

TLT's derivation theorem (0]) is based on the inadequate approach a). In TLT, the 
derivative f'{t) of f{t) is presupposed to exist in the ordinary mathematical sense and the 
L-transform of f'{t) is expressed by (PQ). The theorem for the first derivative then emerges 
from integration by parts: 

L{f{t)} = / f{t)e-^Ut 
Jo 

= sL{f{t)}-fiOl (40) 

Although this kind of mathematical reasoning is formally correct, the result is not LT- 
consistent. The theorem holds only for d-functions f(t) = fd(t), i.e., functions that at 
t = are derivable in the ordinary sense; d-functions are non-LT-consistent because 

u{t)u{t) + uit); cf. m- 

Below, the LT-consistent theorems for derivation and integration are obtained by the 
approach b). These theorems are termed the primary derivation and integration theorem, 
respectively. They hold for derivatives/integrals of the form [u{t) f {t)]^'^'^\ There also 
exist secondary theorems; these hold for /"^^"^(t). TLT's derivation theorem turns out 
to be of the secondary type. (In TLT the distinction between /'•^"•'(t) and [u{t) f (t)]^"^^^ 
is ignored.) The primary theorems for derivation/integration have the same form such 
that they can be unified. The unified theorem can be generalized for non-integer order of 
derivation/integration. The secondary theorems turn out to be virtually irrelevant. 



3.1 The primary derivation theorem 



As the inverse transform ip{t) is a causal and non-ordinary function, its derivatives are 
non-ordinary and causal, as well. By contrast, the formula (0) can within the integral be 
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derived for t in the ordinary sense, and the order of derivation is unhmited, as only the 
function exp(st) needs to be derived. For the first derivative of (f{t) one obtains 

^'(t) = L~'{sFis)} = L~'{sL{f{t)}} = L~'{sL{uit)fit)}}. (41) 

Observing the concatenation theorem (jHT|) . one obtains for the second derivative 

^"(t) = L~'{sLW'it)}} = L-'{s'LWit)}} = L-'{s'L{fm, (42) 

Continuing with this kind of reasoning to obtain higher-order derivatives, and observing 
(j30|l . one obtains the 

primary derivation theorem 

L{[u{t)f{t)p} = s^L{f{t)y, n = 0,l,... (43) 

Equation holds for any type of real function f{t) and for any n for which the inverse 
transform of s"L{/(t)} exists. 

For causal functions fc{t) the primary derivation theorem assumes the form 

L{ft\t)} = s-L{m}; 72 = 0,1,... (44) 

The form ()44|) applies, in particular, to the impulse functions fc{t) = 5^"'\t), and therefore 
Eqs. (jlll are consistent with (j44j) . In contrast to a widespread misconception, the 
reason why the form (j44p holds for impulse functions is not because these are distributions 
but because they are causal. 

The t-domain implications of the operation s"'L{f{t)} are quite different for causal 
functions versus d-functions. As for causal functions fdt) there holds (p{t) = fc{t), there 
follows from the above deduction of (j43|) 

L-'{s-L{m}} = f!:-\t). (45) 

For d-functions the corresponding relationship is considerably more complicated. To 
demonstrate this, functions of the type 

/ud(t) = u{t)U{t) (46) 

are taken into consideration, where fd{t) denotes a d-function. The n-th derivative of the 
so-called ud-function fud(t) is depicted by the 

dud-theorem 

n-l 

fi^it) = Ht)mP = u{t)f!r\t) + E ft'"'\0) ■ 5^''\t)- n = 1, 2, . . . (47) 

v=0 
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The dud-theorem (derivation of ud-function) expresses the n-th non-ordinary derivative 
of u(t)fd(t) by the ordinary derivatives of fd(t), i.e., at the expense of getting impulse 
functions involved as depicted by (jUj). The dud-theorem is explained in Sect. I A. 41 

By (jlZj), utilizing pHjl . the effect of the operation s"L{/d(t)} gets depicted by 

n-l 

L-'{s-L{m}} = n(t)/i")(t) + E ft'~'\0) ■ n = 1, 2, . . . (48) 

As an example, consider the operation sL{cosujt}. From ()48p one obtains 

L'^{sL{cosujt}} = ~u{t) -ujsmujt + dit). (49) 
The same result is obtained from the LT-correspondences ()1331 1134]) . i.e., 

2 

L-^{sL{cosujt}} = L-^\s ■ , ^ A = L-^\l - A = 5 ( t) - u(t) ■ u sin cut. (50) 

The most important domain of application of the primary derivation theorem is de- 
termination of the evoked solution of the inhomogeneous linear DE; cf. Sect. 14.11 

3.2 The secondary derivation theorem 

Although d- functions and their derivatives are non-LT-consistent, the question is legiti- 
mate how the L-transform of fj(^\t) can be expressed by the L-transform of fd{t)- The 
answer is already implied in PH|). From that expression one obtains by L-transformation 
the 

secondary derivation theorem 

L{f!r\t)} = s-Lium - y: ft'-''\o) ■ ^ = i, 2, . . . (51) 

The secondary derivation theorem ()51|) turns out to be identical to TLT's deriva- 
tion theorem. This is a consequence of TLT's original endeavour to provide an operator 
calculus essentially for d-functions; cf. Eq. (jlU)). 

The secondary derivation theorm is not a self-contained derivation theorem as it just 
emerges from application of the primary derivation theorem to d-functions. 

3.3 The primary integration theorem 

According to (j21) the inverse transform's integral function can be depicted by 

^(-i)(t) = L-'{s-'Fis)} = L-'{s-'L{fm = L-As-'L{u{t)fm- (52) 
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As ip(t) is causal and LT-consistent, the integral function Lp^ ^\t) is causal and LT- 
consistent, as well. Thus for the second-order integral there holds 

= L-'{s-'LW^~'\t)}} = L-'{s-'{LWit)}} = L-'{s~'L{fm; (53) 

and so on for higher-order integrals. Utilizing (jHUj) one obtains the 

primary integration theorem 

L{Kt)/(t)](-")} = .-"L{/(t)}; n = 0,l,... (54) 

which for causal functions f(t) = fc{t) reads 

L{ft''\t)} = s--L{Mt)}. (55) 

The operation s~"L{/(t)} is equivalent to n-fold integration of the causal function 
[u{t)f{t)], which implies (n — l)-fold iteration of the t-domain operation 

[n(t)/(t)](-i) = r n(r)/(r) dr = u{t) f f^r) dr. (56) 

J-oo JO 

As u{t)f{t) and the integral functions are causal, the ra-fold integral can be expressed by 
the formula 

[«(t)/(t)](-") = f\t - rr-'fir) dr; n = 1, 2, . . . (57) 

[n — ly. Jo 

The consistency of this formula with the primary integration theorem (jMj) can be verified 
by L-transformation and application of the convolution theorem utilizing the LT- 
correspondence ()130|) . 

The t-domain implications of the operation s^'^L{f{t)} are just as different for causal 
functions versus d-functions as was found for derivation. For causal functions fc{t) one 
obtains 

L-'{s-L{m}} = ft''\t)- (58) 

When f{t) = /d(t) is a d-function such that /d(i) is the n-th ordinary derivative of 
/j""^(t), there holds 

[U{t)f,{t)t'^ = f u{t)U{T)dT = u{t) f U{T)dT 
J-oo Jo 

= um!r"\t) - ft'\0)]. (59) 
By iteration of (j59|) one obtains the so-called 
iud-theorem (integration of ud-function) 

n— 1 Jiy 

/ir^t) = mUit)]^-^ = u{t)ft''\t) - u{t) Y: ft''^'\0) n = 1, 2, . . . (60) 

iy=0 
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For the effect of the operation s "L{/d(t)} one obtains from (jMj) and (lUUj) 

L-'{s--L{m}} = u{t)fi-\t) - u{t) E /1-"^^H0) ■ ^; n = 1,2, . . . (61) 

As an example, consider the operation s'^ L{exp{—at)} . From (j6T| one obtains 

L-i{s-iL{e-»*}} = ■ e-»* + (62) 

a a 

Using the LT-correspondences ()13H 1132]) one obtains the same result, i.e., 

L-'{s-'L{e-'^'}} = L~U- ■ = ^(1 - e-'^*). (63) 

s + a 

3.4 The secondary integration theorem 

By analogy to the secondary derivation theorem, Eq. ()6H) enables for d-functions ex- 
pression of the L-transform of the n-th order integral in terms of the L-transform of the 
d-function itself. By L-transformation of (|6ip. utilizing (|13U|) . there emerges the 

secondary integration theorem 

n-l 

Lift'^Ht)} = s-'-LiUit)} + Y: ft''^'\0) ■ s-^~'; n = 1, 2, . . . (64) 

i/=0 

The secondary integration theorem is complementary to the secondary derivation theorem. 
For n = 1 the two theorems are equivalent; indeed, from ()64j] one obtains 

L{ft'\t)} = s-'L{Ut)} + ft\0) ■ (65) 

Taking into account that by definition fd{t) is the first ordinary derivative of fj^ ^\t), 
Eq. is identical to the secondary derivation theorem and thus to TLT's derivation 
theorem. 

As TLT's derivation theorem (01 En|) is identical to the secondary derivation theorem 
()51|). one would expect TLT's integration theorem to comply with the secondary integra- 
tion theorem ()64|) . However, TLT's integration theorem (0) actually is identical to the 
primary integration theorem (j54j) . The latter theorem holds only for causal functions; cf. 
Eq. (|55|): its application to d-functions - such as is customary in TLT - will in general 
yield erroneous results. 

Thus, it turns out that the alleged mathematical consistency of TLT's theorems for 
derivation/integration is delusive. TLT's derivation theorem holds only for d-functions 
whereas TLT's integration theorem holds only for causal functions. In TLT both theorems 
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are ordinarily used in an untenable way: TLT's derivation theorem is regarded as the 
derivation theorem of LT although it holds only for the non-LT-consistent d-functions. 
TLT's integration theorem is ordinarily utilized for d-functions although it does not apply 
to this type of function. Eventually, these observations explain why there is a formal 
conflict between TLT's theorems for derivation and integration (jUEl). 

The secondary integration theorem is not a self-contained integration theorem as it 
just emerges from application of the primary integration theorem to d-functions. 

3.5 The generalized derivation/integration theorem 

The LT-consistent theorems for derivation and integration, i.e., (j43j) and (j54|) . can obvi- 
ously be unified into one formula, i.e., 

L{[u{t)f{t)p} = s"L{f{t)} n = 0, ±1, ±2, . . . (66) 

Thus, ultimately t-domain integration is in the L-domain merely "reciprocal" to t-domain 
derivation. The t-domain implications of the operation s"L{/(t)} that are outlined in 
Sects. 13.11 13.31 must be observed. 

As a consequence of the formal equivalence of derivation and integration in the L- 
domain description, the theorem can be generalized for non-integer order of deriva- 
tion/integration. For any real r > there holds (Riemann, Liouville, Cauchy) 

D-'Mt)m} = ^ J\t - ry-'m dr; r > 0. (67) 

The operator D denotes generalized derivation/integration, while the negative exponent 
— r indicates r-th-order integration. As the integral in (|67|) is equivalent to the convolution 
f"^ * f{t), Eq. ()67|) can by the convolution theorem (jH^ be L-transformed, and when the 
LT-correspondence ()141|) is utilized one obtains 

L{D~^{u{t)fm = ^Mt-'}-L{f{t)} (68) 

r(r) 

= s-'-L{/(t)}; r > 0. (69) 

Equation (jU^ depicts the generalized integration theorem. This theorem suffices to 
obtain the same kind of generalization for derivation, namely, by concatenation of (j69p and 
(|43p . The sequence of r-th order integration and n-th order derivation yields {n — r)th 
order derivation or integration, depending on whether n > r or n < r JU]. Denoting 
n — r = a one obtains the 

gdi-theorem (generalized derivation/integration theorem) 



L{D"{u{t)f{t)}} = s"L{f{t)y, ae^. 



(70) 
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In the form ()7U|) the theorem holds for any type of real function f{t). For integer values 
of a Eq. (fTOj) is equivalent to (jMj) . 

The concept of fractional calculus, i.e., utilization of derivatives of non- integer order, 
is essentially based on the formula (jEZJ, and this formula holds only for causal functions. 
In the realm of TLT derivatives of non-integer order can not be consistently expressed 
at all, as TLT's derivation theorem neither holds for causal functions nor for non-integer 
order of derivation. When TLT's derivation theorem is used anyway, initial values become 
involved which merely are a nuisance. It is by the present approach - which provides for 
the generalized theorem - that Laplace transformation becomes a consistent and 
efficient tool for doing fractional calculus. 

Examples for the application of (f7n|) . i.e., for a = 1/2, are listed in Sect. IA.5I 
3.6 Summary 

LT-consistent theorems for t-domain derivation and integration are obtained by deduction 
from the behavior of the inverse L-transform. These theorems - the primary theorems 
- are formally congruent and thus can be unified. The unified primary theorem is gen- 
eralized for non- integer order of derivation/integration. The generalized theorem (the 
gdi-theorem) accounts in the L-domain for both integer-order and non-integer-order of 
t-domain derivation/integration. As the inverse transform is causal the primary theorems 
and the gdi-theorem hold for causal functions. 

The respective secondary theorems emerge from application of the primary ones to 
the L-transforms of ordinary derivable functions (d-functions). TLT's derivation theorem 
is identical to the secondary derivation theorem. This theorem can not be regarded as 
the derivation theorem of Laplace transformation because it applies only to the class of 
d-functions, which is non- LT-consistent. TLT's integration theorem is identical to the 
primary integration theorem and therefore holds only for causal functions; its utilization 
for d-functions - such as in TLT - is untenable. 



4 Solution of linear differential equations 

Below, the new approach to the solution of linear DEs by LT is demonstrated for a fairly 
general form of the ordinary linear DE, namely, 

N M 

$:an2/(")(t)= Y.braX^'^Hty, iV = 1, 2 . . . ; M = 0, 1 . . . (71) 

n=0 m=0 

The coefficients presupposed to be real constants. Two functions are involved, 

namely, the excitation function x{t) and the response function y{t). The form (f7T|) is 
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more general than that accounted for by the TLT method, as in ()71|) derivatives of the 
excitation function are allowed to be included, i.e., corresponding to M > 0. 

A clear distinction is made between the evoked solution and the spontaneous solution of 
the inhomogeneous DE |13^. The evoked solution (the pertinent system's evoked response 
to x{t)) is that part of the total response that is elicited by the excitation function x{t) 
alone. The spontaneous solution (the system's "spontaneous" response), if it exists, is 
ascribed to a pre-excited initial state of the system. Mathematically, the evoked solution is 
the DE's total (general) solution minus the general solution of the pertinent homogeneous 
DE. The spontaneous solution is identical and synonymous to the general solution of the 
homogeneous DE. 



4.1 The evoked response: 

Particular solution of the inhomogeneous DE 

For the excitation function x{t) to be LT-consistent it must be defined as a causal function. 
As a consequence, its derivatives have the form [u{t)x{t)]^"^\ Another consequence is that 
the evoked response ye(t) is causal as well; thus, its derivatives implicitly assume the form 
[u{t)yc{t)]^'^\ Therefore, conversion of ()71|) into the L-domain is governed by the primary 
derivation theorem ()43|1 . The L-domain representation of (ffT| reads 

N M 

J2 anS"L{y,(t)} = J2 b,,,s^^L{x{t)}; AT = 1, 2 . . . ; M = 0, 1 . . . , (72) 

n=0 m=0 

and one eventually obtains the evoked solution 

yeit) = L-'\^^^-L{x{t)}]. (73) 

In contrast to the TLT method, the evoked response is by (fTSj) obtained without 
involvement of initial values, i.e., additional constants. To obtain the evoked response 
one does not have to pretend that the pertinent system is in a non-preexcited initial 
state. With respect to the evoked response the system's initial state is irrelevant. 

For the operation L~^, i.e., inverse L-transformation, there exist two alternatives. The 
first of them is taking ()73|) at face value, which implies that the L-transform of x{t) has 
to be included. The second alternative exploits the convolution theorem ((221), whereby 
inverse L-transformation can be confined to the first factor in (fTSj) : 

hit)=L-'\^^]. (74) 

The function h(t) depicts the evoked response to the excitation function x{t) = S{t), i.e., 
the impulse response. As (|73j) has the form 

y,{t)=L-'{L{h{t)}-L{x{t)}}, (75) 
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by ()32|) there emerges the famihar convolution formula 

ye{t) = u{t) f hit - t)x{t) dr. (76) 



For example, for = 1, M = one obtains the impulse response 

h{t) = L-A I = M(t)/?e-°*; a = ao/ai; /5 = 6o/ai. (77) 

{ao + ais J 

Utilizing (f77j) and (fTBj) . one obtains the 

evoked solution of the first-order linear DE 

y^it) = u{t)P ra;(r)e-"(*-") dr = n(t)/3e-"* f x{T)e''^ dr. (78) 
Jo Jo 

Though this solution is entirely expressed in the t-domain, it nevertheless is based on LT, 
which implies that ()78|) is correct only for causal excitation functions x{t). 

It may be noted that it is the independent LT-based expression of the evoked solution 
(I73p that provides a solid basis to operator calculus as it is customarily employed in the 
theories of linear control systems and electrical circuits. 



4.2 The spontaneous response: 

General solution of the homogeneous DE 

When the evoked response is by LT determined as just described, for the solution of the 
homogeneous DE still any method available can be chosen. From the present approach 
there emerges a new LT-based method which exploits the facts that 

a) the spontaneous response ys{t) and its derivatives are d-functions; their L-transforms 
have the form L{yj^''\t)} = L{u(t)?/(")(t)}; 

b) the homogeneous DE can by the dud-theorem be converted into an equivalent inho- 
mogeneous DE. 

To get the homogeneous DE made up for conversion into an equivalent inhomogeneous 
DE, it is multiplied by u{t); this yields 

aou{t)y,{t) + aiu{t)y'^{t) + ... + aNu{t)yf\t) = 0. (79) 

Because of LT's ambivalence, the L-transform of (|79|) is identical to the L-transform of 
the original homogeneous DE. Application of the dud-theorem (I47|) converts (|79|) into the 
form 

aou{t)ys{t) + ai{[u{t)yM' - VsW ■ S{t)) + a2{[u{t)y,{t)]" - y'M ■ S{t) - y^ " S'{t)) 

N~l 

+ ... + a;v([«(t)l/s(t)](^) - E yf-'-'\0) ■ 5(^)(t)) = 0. (80) 

u=0 
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Equation (jHUj) is equivalent to (f?^ . However, (IHUj) actually is an inhomogeneous DE. The 
impulse functions included in ()80p play the role of virtual excitation functions jT^j. This 
becomes particularly apparent when (j80|l is expressed in the form 

N N-l 

Y.anHt)yMP=T.^,^^'\t)^ (81) 

n=0 /J=0 



where the coefficients are depicted by 



Af-l- 



E 

u=0 



a 



(82) 



From (jHT|) the solution of the homogeneous DE can be obtained as an evoked solution, 
i.e., as described in Sect. 14.11 One eventually obtains 

»(%.(*) = i- 1 I . (83) 



The notation u{t)ys{t) is not redundant, because Usit) is a d-function. Equation (j83|) 
in fact depicts the ewoA;ec? response of a virtual system, i.e., to the excitation function 6{t). 
The transmission function of that virtual system is defined by the quotient in (jHH|) . whose 
numerator essentially originates from the impulse functions contained in (jHT|) . 

Finally, the spontaneous solution ys{t) itself can by extrapolation into t < be ob- 
tained from (jSH|) . As according to the causahty theorem (jHU)) the operation L^^ on the 
right side of (|83p yields either explicitly or implicitly a function of the form u{t)ys{t), that 
kind of extrapolation is equivalent to cancellation of u{t) on both sides of ()83|) . i.e., after 
inverse L-transformation. 

For example, for = 1 one obtains from (j83p and (j82p 

u{t)y,{t) = L-'l^^yf^]=u{t)yMe-'''; a = ao/a,. (84) 
[ ao + ais J 

By cancellation of u{t) there emerges from ()84|) the 
spontaneous solution of the first-order linear DE 

ysit) = i/s(0)e-"* for - oo < t. (85) 



The described deduction of fjHHj) ultimately constitutes a new method for obtaining 
the general solution of the linear homogeneous DE of finite order N, i.e., by Laplace 
transformation. The method involves a) determination by ()82|) of the coefficients c^, i.e., 
from the DE's coefficients a„; b) insertion of the coefficients into (j83|) : c) inverse L- 
transformation; and d) cancellation of the factor u{t). The arbitrary constants which 
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invariably get involved are provided by the initial values at t = of the solution ys{t) itself 
and of the latter's — 1 derivatives. The method is straightforward; yet, mathematical 
intricacies may occur in step c), i.e., inverse L-transformation. 

As a pragmatic alternative to the described deduction, the formula (fHH|l may be ob- 
tained in a more immediate way, namely, by using for LT-conversion of the homogeneous 
DE the secondary derivation theorem ()51|). This option is enabled by the fact that the 
secondary derivation theorem accounts for the combination of the dud-theorem with the 
primary derivation theorem; cf. Sect. 13.21 As the secondary derivation theorem is identical 
to TLT's derivation theorem, this option explains why TLT's total solution of the inho- 
mogeneous DE is quite similar to the new total solution; see the following two sections. 

4.3 The total solution 

Assuming that for the solution of the homogeneous DE the method just described is 
chosen, the total solution of (f7T|) can be compactly depicted by a formula, i.e.. 

The first term on the right side of (|H^ depicts the evoked response and is identical to 
(fTSj) . The second term depicts the spontaneous response by one single expression; this 
term is equivalent to the combination of Eqs.()83|) and ()82|). The operator indicates 
inverse transformation by (jSJ. The operator Lj^ indicates inverse transformation by 
followed by extrapolation into t < 0, i.e., cancellation of the factor u{t). 

When the inverse L-transform is looked up from a TLT-based table of LT-correspon- 
dences it must be observed that in those tables the t-domain functions are depicted in 
the non-causal form. In Sect. I A. 71 examples are listed of the LT-consistent notation of 
t-domain functions. 

It must be kept in mind that LT-based solutions of linear DEs hold only for causal 
excitation functions x{t) = u{t)x{t). The solution for the "steady-state" may be obtained 
by asymptotical approach, i.e., for t oo. 

An example for application of ()86p is depicted in Sect. IA.61 

Finally, it should be noticed that for = 1, M = the LT-based solution of (f7T|) can 
be entirely expressed in the t-domain. By superposition of Eqs. (f75j) and (jH3j) one obtains 
the 

total solution of the first-order linear DE 

y{t) = y,{t) + y,{t) = u{t)l3e-''' f a;(r)e"" dr + ?/s(0)e-"*; a = ao/av, j3 = bo/ai. (87) 

Jo 
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Thus, for the special case that in (jTip there is = 1, M = one does not actually 
have to do L-transformation at all. The solution is reduced to evaluation of the integral 
contained in ()87|1 and it holds for causal excitation functions x{t) = u{t)x{t). 



4.4 TLT's initial- value conflict 

One of the most annoying deficiencies of TLT is potential interference of the DE's evoked 
solution with the spontaneous solution. The danger of this kind of interference to occur 
emerges from the fact that the DE's solution obtained by TLT includes an inappropriate 
type of initial values. This is why the phenomenon is termed the initial- value conflict. 

When the general solution of a linear DE of the form (j71|) is worked out by means of 
TLT - which is possible for M = - one eventually obtains the formula 

This formula is quite similar to (jSUj) . i.e., for M = 0. (The constant Bq is in fjSH|) included 
just for formal compatibility with (j86p : one may in both equations assume Bq = 1.) Letting 
aside the difference in inverse L-transformation that is apparent by comparison of (jHHjl 
to (jHUj) . there remains the crucial difference that (jHUj) correctly contains the initial values 
of the spontaneous solution, ys'^\0), whereas (j88|l contains the initial values of the total 
solution, y^^^O). The latter, inadequate kind of initial values inevitably emerge from 
LT-conversion of the DE ()7H1 by TLT's derivation theorem. 

The potential initial-value conflict arises from the fact that the evoked and the spon- 
taneous solutions are superimposed in y{t), such that 

^H(0)=yf(0)+yf (0); ^. = 0,1,... (89) 

The initial values of the evoked solution, y^'^\0), depend on x{t) and on the DE's co- 
efficients Cljj, clS IS depicted by (f7H|) . If there happens to be y^'^\0) = 0, i.e., for all 
z/ = 0, 1, . . . , - 1, then there holds y'^''\0) = yi^^O), and (jHHD will yield the same re- 
sult as ()86p. i.e., except for the difference in inverse L-transformation, and for M = 0. 
However, for many types of system and of excitation function yi'^\0) may turn out to 
be different from null, i.e., for at least one value of u. If this occurs then the result ob- 
tained by TLT is incorrect; the initial values y^'^\0) can no longer be freely chosen, i.e., 
independent of the excitation function; cf. Sect. lA.ll 



Obviously, (|88|) can easily be "patched", i.e., by arbitrarily substituting the values 
y^'^\0) for y^'^\0). Most remarkably, there exists another, less trivial kind of patch, namely, 
arbitrary substitution in (|88p of the values y^'^\—0) for y^'^\0). The reason why the 
latter kind of patch works around the initial-value conflict originates from the fact that 
the evoked response and its derivatives are causal (cf. Sect. 14. such that there holds 
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yi'^\—0) = 0, i.e., even when y^'^\0) ^ 0. As the spontaneous response and its derivatives 
are d-functions, there holds y^'^\—0) = y^'^\0). Thus, one obtains from (jH^ 

l/(^)(-0) = yi'^O) for z/ = 0, 1, . . . , iV - 1. (90) 

The identity ()90p holds for any type of linear DE and of excitation function. From this 
identity there follows that substitution in (j88|) of the initial point t = —0 for t = 
renders (jHH|) compatible with (|HH|l . such that TLT's solution becomes essentially correct, 
i.e., except for the difference in inverse L-transformation. The "empirical" finding that 
by this kind of patch TLT's initial-value conflict actually is worked around, becomes in 
this way explained. 

It should be noticed that the latter kind of patch is just as arbitrary as the former. 
Even if one had no doubt at all about TLT's consistency, the mere fact that TLT's solution 
fl88|) in general needs to be patched should have elicited that kind of doubt. One can not 
in earnest be satisfied with either patch's working around TLT's initial-value conflict. 

4.5 Summary 

The total solution of the inhomogeneous linear DE is obtained by two independent algo- 
rithms, i.e., one for the evoked response, another for the spontaneous response. As both 
the excitation function and the evoked response function are causal, the derivatives of 
these functions assume the form [u{t)x(t)]'^"^^ and [u{t)yc{t)Y'^\ respectively. When the 
inhomogeneous DE is by the primary derivation theorem converted into the L-domain one 
obtains a linear algebraic equation that does not contain any arbitrary constants. This 
equation is solved in the familiar way to obtain the evoked solution yeit) which is defined 
for — oo < t. 

For the determination of the spontaneous response ys{t), i.e., solution of the pertinent 
homogeneous DE, the fact that the response function and its derivatives are d-functions 
is a priori taken into account. The L-transforms of these functions assume the form 
L{y^'^\t)} = L{u{t)y^'^\t)} . As a consequence, the homogeneous DE can be multiplied 
by u{t) without affecting its L-transform. The homogeneous DE such modified can by 
the dud-theorem be converted into an equivalent inhomogeneous DE. The L-transform of 
the latter DE depicts a virtual linear system whose impulse response equals for t > the 
spontaneous solution of the original homogeneous DE. By extrapolation into t < one 
obtains from that impulse response the spontaneous solution y^it) for — oo < t. 

For M = 0, i.e., if the excitation function does not include any derivative, the total 
solution obtained by TLT is formally similar to the correct solution depicted by the new 
method. However, the TLT-based solution suffers from an initial- value conflict. This 
conflict can be worked around by manipulation of the initial values. The necessity of 
patching TLT's solution is another symptom of TLT's deficiency. 
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5 Concluding remarks 

From the theorems pointed out in Sect. |21 it should be apparent that the key to the 
consistent LT-theory hes in the hidden bilaterahty of the transformation integral (^. By 
strictly taking that bilaterality into account there emerges a structure of LT-theory that 
is fundamentally different from that of TLT. Yet, many of the new theory's theorems and 
expressions are familiar from TLT. Pronounced formal differences from TLT's expressions 
occur, e.g., in the context of the theorems for derivation and integration. The difference 
between the new formula ()86|) for the general solution of the linear inhomogeneous DE 
from the corresponding formula ()88p obtained by TLT, though at first sight marginal, 
actually is crucial, namely, for resolution of TLT's inherent initial-value conflict. 

Mathematical consistency of the new approach to LT is achieved at the expense of non- 
ordinary t-domain functions being crucially involved - primarily, the inverse L-transform 
(p{t), which includes the connect function Uo(t) and, possibly, its (non-ordinary) deriva- 
tives. This notion brings to mind the aspect of mathematical rigour. If one maintains that 
there is no rigorous mathematical account for non-ordinary functions and non-ordinary 
derivatives then the inverse L-transform can not in general be rigorously accounted for 
and a rigorous theory of Laplace transformation can not exist at all. 

Rigorous treatment of non-ordinary functions may be achieved by invoking the theory 
of distributions. However, what makes LT-theory consistent in the first place is adequate 
incorporation of non-ordinary functions and non-ordinary derivatives. Mere attachment 
of distribution theory to TLT can not cure TLT's inconsistency. Mathematically rigor- 
ous treatment of non-ordinary functions and derivatives may be achieved separately and 
subsequently - which, historically, actually has happened to the delta-impulse. 

Making in this sense a distinction between mathematical consistency and rigour, one 
can say that the new approach to Laplace transformation outlined in the present article 
is consistent, even though it may not fully qualify as mathematically rigorous. This 
achievement is distinctly preferable to inconsistency disguised by fake mathematical rigour 
- which virtually is what is offered by Doetsch jH |S] and many others. 

Once one has developed some scepticism about TLT one may be prepared to realiz- 
ing that Doetsch's books are written in a somewhat defensive and dogmatic style, and 
that Doetsch devoted considerable portions of text to reasoning by plausibility and to 
explaining away apparent discrepancies. From these observations it may be concluded 
that Doetsch virtually was aware of TLT's inherent problems. Indeed, in the preface to 
[2j he suggested the need for a fundamental redesign of TLT. However, until the end of 
his life (1977) he stuck to the original layout of TLT, i.e., with some distribution theory 
as a complement. In the two volumes of his Handbook OEl, Doetsch assembled a wealth 
of mathematics on LT - which, however, may distract from TLT's fundamental flaws. 

To obtain a more thorough understanding of TLT's unfortunate history it is probably 
helpful to take notice of Doetsch's biography jHUHI- 
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A APPENDIX 



A.l Failure of TLT: An example DE 

TLT promises to provide a method for obtaining the general solution of the linear in- 
homogeneous DE, i.e., for any kind of excitation function x{t). Below, this promise is 
challenged by an example, i.e., for the first-order DE 

ay(t) + y'{t) = x sin ojt, (91) 

where xsinujt denotes the excitation function and y{t) denotes the response function. 
From the theory of linear DE's one obtains the general solution 

X 

y{t) = -^-^ — -{asinut - u cosut) + ys{0)e~"-\ (92) 

where ys{0) is an arbitrary real constant which specifies any particular solution of the 
pertinent homogeneous DE. 

By the TLT method one obtains the solution 



^s^ + tu^ s + aj [s + aj 

^ 'ooe'"^ + asinut - ujcosut) +y{0)e~''^; for t > 0. (93) 



The result fl93|) differs from the correct solution ()92|) in two significant respects: 

a) As compared to the correct evoked response 

X 

yJt) = — -(asinut — u cosut) (94) 

-|- 

Eq. ()93|) includes an extra decaying exponential. 

b) As compared to the correct spontaneous response 

ys{t) = i/s(0)e-'^* (95) 

Eq. includes the initial value y{0) of the tota/ response y{t) = ydf) +ys{t) instead of 
that of the spontaneous response alone, i.e., |/s(0). 

The solutions ^ and ^ both are consistent with the DE dHD). Yet only ^ 
represents the DE's general solution. This becomes apparent when one attempts formal 
reconciliation of (j93j) with (j92p . Reconciliation can only be obtained by suitable choice of 
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a particular "initial state". For instance, for ys{0) = Eqs. can only be reconciled 

by setting 



Hence, y{0) is not really a free parameter of the spontaneous response. While the consis- 
tency of with suggests that the TLT-based solution is correct, comparison with 
()92|) reveals that actually is only a particular solution of (j^Tj) and that y{0) is not a 
freely assignable constant. 

In point of fact, the correct solution of (|9ip can by L-transformation not be obtained at 
all, because the L-transform of x sin cut is just an alias for L{u{t) ■ x smut}; cf. Sect. 12.41 
The evoked part of the TLT-based solution (j^J^ inevitably depicts for t > the response 
to u{t) ■ xsinut instead of to xsincot. The above observations reveal that LT actually is 
sensitive to the difference between x(t) and u{t)x{t), i.e., if there is x(t) ^ u(t)x{t). 

From a more general point of view, it should be appreciated that any "causal" physical 
system which is governed by a linear DE of finite order has a "memory" which becomes 
manifest in the length of the system's impulse response. This notion suffices to explain 
that the evoked response ye{t) to an excitation function x{t) in general depends on the 
behavior of x{t) for —oo < t fUE]. When the response is determined by LT, this fact 
is automatically accounted for, namely, by the implicit bilaterality of (H)); cf. Sects. 12.11 
12. 4[ IA.21 The bilaterality of ip(t) is indispensable for the evoked response to be correct, 
i.e., whether or not (f{t) is causal. L-transformation by (P) just provides for (f{t) being 
always causal. As a consequence, the evoked response complies with x{t) only if x{t) was 
a priori defined to be causal. 

A. 2 Analysis interval versus integration interval 

In Sect. 12.11 it is pointed out that the t-domain interval which actually is encompassed 
both by Fourier- and Laplace-transformation invariably is infinite, i.e., — oo < t < +oo. 
This interval is termed the respective transformation's analysis interval; it has to be 
distinguished from the integration interval, i.e., the interval of t which is determined by 
the transformation-integral's limits, e.g. in Eqs. and (jH)). The relationship between 
these two kinds of interval can be made apparent by regressing to discrete Fourier- and 
Laplace-transformation. The discrete transformations can be deduced from the Fourier- 
series representation of a real function f{t), i.e.. 



oo 



oo 



/(t)=ao+E 



On COS Unt + X] 
n=l 



sincUnt; n integer; 



(97) 



n=l 



where 




(98) 



Ernst Terhardt: Laplace transformation updated 



29 



2 /•+^/2 

On = — / f{t) COS dt] (99) 

1 J-T/2 

6n = -/ /(t)smcunt dt; (100) 
2vrn 

= (101) 

The function f(t) such represented is periodic with the period length T = l-ajuox. 

This set of formulas can be regarded and employed as a complementary pair of trans- 
formations ^T]. By Eqs. ()98II100|) a particular section of f{t) becomes transformed into 
an infinite set of real coefficients {a^, b^} each of which pertains to a discrete frequency u^^. 
The section of f{t) that extends from — T/2 to +T/2 plays the role of an analysis interval, 
T. The function f{t) needs to be defined for —T/2 < t < +T/2. Inverse transformation 
is achieved by (jHIj). In general, the inverse transform matches the original function f{t) 
only for —T/2 < t < +T/2. By the notation chosen in Eqs. ()98lll()0j) the center of the 
analysis interval becomes implicitly denoted t = 0. 

To elucidate the implications of this approach it is helpful to express the above formulas 
in complex notation. Using Euler's formula, one obtains from Eqs. (|99( llOOp 

an ± i6n = - / /(t)e±'"'^* dt; n = 1, 2, . . . (102) 

1 J-T/2 

By definition of the discrete Fourier transform 

T(a;n,T) = |(an-i6„) (103) 



one obtains from ()102|) 



r+T/2 

F{uj^,T)= I /(t)e-^^"* dt. (104) 

T/2 



Equation (jHZj) can be expressed in the form 

oo oo 

fit) = ao + i E(«n - i&n)e^""* + I E(«n + i&n)e-'""*. (105) 

n=l n=l 

Observing that by (jnHMl) there holds Oq = T(0, T)/T, one obtains from (ITOHD and ^TU^ 
for n = ...,— 1,0,+1,... the formula for inverse transformation 

1 +00 

^(t) = 7^ E ^(^n, T) ■ e'-"* = fit) for - T/2 < t < +T/2. (106) 

n=— oo 



The transformation ()104|) is governed by, and confined to, the analysis interval T. This 
interval can either be regarded as preset or as being determined by the choice of ui and 
thus of the spacing of the discrete frequencies Un, cf. fll()l|l . As the center of the analysis 



Ernst Terhardt: Laplace transformation updated 



30 



interval is for convenience denoted t = there is only one way in which the limits of 
the integral ()104j) may become different from ±T/2, namely, by f{t) being null within a 
sub-interval of T that borders t = — T/2 and/or t = +T/2. Vice versa, when one or both 
of the integral's limits are different from ±T/2 this invariably indicates that the function 
f{t) contains a null interval. In particular, when f{t) is replaced with the causal function 
fc{t) = u{t)f{t) the discrete Fourier transform gets expressed by 

f+T/2 /■T/2 /■T/2 

FK,r) = / M(t)/(t)e-'-"* dt = / /(t)e-^-"* dt = / /c(t)e-^--* dt. (107) 

J-T/2 JO JO 

The difference from —T/2 of the integral's low limit neither affects the analysis interval 
T nor does it make the transformation unilateral; it just indicates that f{t) is causal. 

The discrete variant of Laplace transformation emerges from p()4j) by including in the 
integral the factor exp(— crt) and by setting the integral's low limit to t = 0. This yields 

LT{f{t)}= f{t)e-'-Ut, (108) 
Jo 

where Lt denotes discrete L-transformation, and 

Sn = 0" + iuj^. (109) 

Regarding the relationship between analysis interval and integration interval, the same 
applies as was just pointed out above: The integral's ()108p low limit t = indicates that 
the function actually transformed is causal while the low limit marks the center of the 
analysis interval. 

The inverse discrete L-transform ip{t) is determined by the expression 

1 +00 

^{t) = L^'{Lr{fm} = 7f; E ^T{/(t)}-e-*; -T/2<t<+T/2. (110) 

n=—oo 

When ()108|) is employed for transformation of a non-causal function f{t), the inverse 
transform ip{t) still is defined for —T/2 < t < +T/2; however, validity of (f{t) = f{t) is 
confined to the interval < t < T/2. In this sense, the transformation is for non-causal 
functions inconsistent. For causal functions f{t) = fc{t) the transformation is consistent, 
as there holds (p{t) = /c(t) for -T/2 < t < +T/2. 

Laplace transformation as defined by (ITJI21) emerges from (110811110]) by letting T ^ oo, 
implying that the spacing of analysis frequencies becomes infinitesimally small such that 
the sum ()110|1 becomes an integral, i.e., Q. The transformation integral's (0) low limit, 
i.e., t = 0, marks the center of the infinite analysis interval. 

A. 3 Redundancy of multiplication by u(t) 

When the unit step function is defined without inclusion of the connect function uo{t), 
i.e., u{t) = for t < 0; u{t) = 1 for t > 0, there evidently holds M"(t) = u{t) (n = 1, 2, . . .). 
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Hence, multiplication of u{t) by itself is redundant. As u{t) actually includes 'Uo('^) it must 
be verified that this kind of redundancy also holds when u{t) is defined according to (fT^ . 

For = u{t) to hold it is necessary and sufficient that u^lt) = UQ{t). As was 

outlined in Sect. 12 .31 the unit connect function uo{t) is sufficiently characterized by saying 
that it is an infinite set of real numbers {0 ... 1} that exists at t = and fills the gap 
which otherwise exists between u{—0) = and u{+0) = 1. From this definition it follows 
that any integer power of UQ{t) is characterized by precisely the same criterion: u^it) is 
for any n = 2, 3, . . . also an infinite set of real numbers {0. . . 1}. Thus, any power of Uoit) 
can be renamed UQ{t). Therefore, there actually holds the identity 

tt"(t)=M(t); n = l,2,... (Ill) 

Another important question is whether multiplication by u{t) of a derivative of u{t), 
i.e., of an impulse function, is also redundant such that there holds 

= M(t)(5("-^)(t) = m("H^) = '^^"''Hi); n = l,2,... (112) 

The identity (jll2j) turns out to be an inevitable consequence of the fact that the 
impulse functions are LT-consistent (Sect.|2I21 Eq. (^2))), in combination with the causality 
theorem When in ^ one lets f{t) = 5(")(t) one obtains from ^ and ^ 

L-^{L{6^"^{t)}} = 5(")(t) = M(t)5(")(t). (113) 

Thus, the identity ()112|) actually is implied in LT's fundamental definitions and fea- 
tures. The identity can be made plausible by characterizing 5^^''{t) as a set of pseudo- 
functions the type of which depends on n and which only at t = are different from 0. 
Then it may be concluded that multiplication of 5^''^\t) by Uo{t) does not alter the type 
of function denoted S^'^\t) such that, indeed, u{t)5^^\t) can be renamed S^'^\t). If more 
mathematical rigour is desired, the theory of distributions may be invoked. 

In Sect. 12.41 the identities (|111|) and (|112p are subsumed by the unit-step redundancy 
theorem (j21I). 



A. 4 The dud-theorem 

The dud-theorem depicts the n-th derivative of the ud-function /ud(i':) such as defined by 
PUj). as follows. 

With regard to the definition (|14j) of the unit step function the first derivative can be 
expressed by 

fLit) = Ht)Uit)]' = for t < 0, 

= /d(0)M'o(t) fort = 0, 

= /^(t) fort>0. (114) 
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To obtain the derivatives of higher order, it is helpful that ()114j) can be converted into a 
more convenient form, as follows. Utilizing the definition (fT^ of u{t), one can replace the 
third line on the right side of ()114j) with the expression 

wW/dW-/d(0)«o(t) for -oo<t. (115) 

This converts ()114|1 into a superposition of causal functions which are known to be LT- 
consistent. As these terms are LT-consistent they can be replaced with the pertinent 
inverse L-transforms without affecting the validity of /ud(^)- This has the effect that the 
term /^(O)Mo(t) gets eliminated from pi5|) . As a result, f^^^it) can be expressed in the 
form 

fUt) = <t)m + WXit) for - oo < t. (116) 

By comparison of pi6|) to ()46p the scheme becomes apparent according to which the 
second derivative emerges from the first, the third from the second, and so on. For in- 
stance, when flll6|) is derived to obtain the second derivative of fndit), the term u{t)f^(t) 
gets converted into u{t)f^{t) + /^(O)MQ(t); and the second term /d(0)Mo(t) gets converted 
into /d(0)Mo(t). By this scheme one eventually obtains for the n-th derivative the expres- 
sion 

/id^W = u{t)f!['\t) + ft'\(^K{t) + ft"\0)<{t) + ... + M0W^\ty, n = 1, 2, . . . 

(117) 

Equation (jll7j) is by (jTTj) expressed as the dud-theorem. 
A. 5 Derivatives of order 1/2 

From the derivatives of non-integer order of causal functions those of order 1/2 are of 
particular interest fTU' . Below, a number of derivatives are listed which were determined 



by the gdi-theorem (fTUI) . using the table of LT-correspondences Sect. IA.7I 

D^/^u(t)} = L-^\y^--} = u(t)-^ (118) 

D^^uit) .V-t} = L-'{V-s-^} = u{t) . ^ (120) 

D^/2{n(t) • e*} = L-^{v^ ■ = u{t) ■ [l/vGd + e* ■ erf (v^)] (121) 

D'/'{uit).\nt} = L-Hv^- ^^^^^^1 = u{t) . (122) 
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D'/^{u(t) ■ Jo(2v^)} = L-'\^ ■ e-"/^| 



u{t) 



cos 2 vat 



TTt 



(124) 



A. 6 Linear DE that includes a derivative 
of the excitation function 



When the hnear first-order DE is modified into the form 

ayit) + y'it) = x'{t), 



(125) 



the response y{t) denotes the electrical current through the capacitor of the electrical 
circuit shown in Fig. 3, while x{t) denotes the current exerted on the circuit by the source 
Q. The constant a equals a = 1/{RC). 



or r 



-^i§- 3- Electrical circuit which is accounted for by Eq. (|125j) . x{t): electrical 
current from source Q; rj{t): voltage at condenser C; y{t): current through C 



The DE ()125p is of the type which by Doetsch was banished from his theory because 
it includes a derivative of the excitation function. In terms of (f7T|) this DE corresponds 
to = 1, M = 1. 

The new method, i.e., the formula (j86p . provides for an algorithmic, straightforward 
solution. Assuming, as an example, x{t) = xu{t), and taking into account the parameters 
N = M = 1; aQ = a; ai = 1; Bq = 0; bi = 1, one obtains 

X T _ir?/s(0) 



s + + a- 

[xu{t) + VsiO)] ■ e""* for - oo < t. 



(126) 



Any particular initial state can be freely accounted for by setting ys{0) accordingly. The 
two components of the total solution are illustrated in Fig. 4. 



X- 




rs(0P 







i " 








Fig. 4. The two components of the solution of 
()125|) . Left: evoked response for x{t) = xu{t). 
Right: spontaneous reponse. This solution is de- 
picted by Eq. (tT^ 
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A. 7 LT-correspondences 



Below, LT-correspondences are listed that are used in the present article. The t-domain 
functions are identical to those included in customary tables of TLT, except for the factor 
u{t), which is required to make the correspondences valid in both directions. Notice that 
the first correspondence ()127p holds only in one direction. 

1 1/s (127) 

u{t) ^ 1/s (128) 
5(")(t) =M("+^)(t) =4"+i)(t)^s" (n = 0,l,...) (129) 

uit)-t^^^ (n = 0,l,...) (130) 



u(t) . e-"' ^ ^— (131) 

^ ^ s + a ^ ^ 

^(l_e-*)^_i_ (132) 

a ^ ' s{s + a) ^ ' 

u{t) ■ sin ujt ^ , ^ ^ (133) 



g 

u(t) ■ COS ut^— (134) 

^^^^ „ (a;e~°* + asinut - u cos ut) ^ — (135) 

+u!^ [s + a) [s^ + a;2 j 

u{t) ■ Vt ^ (137) 
2s 



uit) ■ e* ■ erf (Vt) ^ (138) 

(s - l)v^ 



1 

u{t) ■ In t ^ ~ ^ (Ce = 0.577215 . . .) (139) 



s 



u{t) ■ ^ ^ -.[^ ■ (In 4s + Ce) (140) 
V s 

«W ■ ^ ^^^^TTT^ (r G r > -1) (141) 

4£ . e-V(4t) ^ 1 . e-av/i (142) 



u{t) ■ a ■ t-^'^ ■ e-'^'/(4*) ^ 2v^ ■ e''^^^ (143) 
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^ • cos 2 Vai ^ ^ • e-"/^ (144) 
u{t) ■ Jo{2Vai) ^ - ■ e-"/^ (145) 
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